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By Manuel Stein and J. Lyell Sanders, Jr. 


'SUMMARY 


A method is presented for obtaining influence coefficients for thin 
low -aspect -ratio wings of built-up construction. Chordwise deflections 
are assumed to be parabolic (or linear) and the principle of minimum 
potential energy is used in conjunction with difference equivalents to 
obtain an appropriate set of equilibrium equations. Symmetric and anti- 
symmetric load-support conditions are considered. A simple method is 
given for taking account of attachment to a flexible fuselage. The com- 
putations involved are for the most part organized into matrix form so 
as to be suitable for high-speed computing machines. The matrices 
needed at the beginning of the computing process can be set up directly 
from the data of the wing design. 


INTRODUCTION 


The simple beam theory based on the concept of an elastic axis with 
uncoupled bending and twisting distortions has long been successfully 
applied to high-aspect-ratlo unswept wings. However, the assumptions 
upon which the simple beam theory is based invalidate its application 
to low -aspect-ratio wings. In the case of a low -aspect-ratio wing, 
bending and torsion are not uncoupled; among other things, there appears 
the additional complication of chordwise bending. Thus, the development 
of an adequate load-deflection analysis for use in predicting the aero- 
elastic behavior of an airplane with low-aspect-ratio wings has became 
an important structural problem. 

Several approaches to the solution of this problem have appeared in 
the recent literature. In a method introduced by Levy (ref. 1), the wing 
is idealized into a network of crisscrossing beams to represent ribs and 
spars with four-sided torsion boxes attached at their four comers to the 
intersections of the beams to represent the shear -carrying capacity of 
the skin. (See fig. 1.) An effective width of the skin is allotted to 
the beams to account for its direct-stress-carrying capacity. A stiffness 
matrix (or inverse influence-coefficient matrix) is found for the composite 
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structure by properly combining the stiffness matrices of the various 
components. The inf luence -coefficient matrix of the composite structure 
is obtained by inverting the stiffness matrix of the composite structure. 
The method is meant to be applicable to thin wings with thin cover sheets. 

In the "wide beam" theory introduced by Schuerch (ref. 2), the wing 
is idealized into a bundle of alternating simple beams and torsion tubes 
running in the spanwise direction. (See fig. 2.) The beams represent 
the spars and axial-load-carrying capacity of the skin. The torsion tubes 
represent the shear-carrying capacity of the skin. The bundle is tied 
together by ribs which are assumed rigid. The individual beams and tor- 
sion tubes are assumed to behave according to simple beam theory. The 
equations of moment and torque equilibrium for the idealized structure 
are two simultaneous differential equations (both of the fourth order) 
for the deflection and rotation of the cross sections. These equations 
are not uncoupled as in simple beam theory, and several new section prop- 
erties appear. This method, presumably, is meant to be applied to thin 
wings with thin cover sheets . In reference 3 Schuerch and Freelin have 
presented another method which is somewhat like that of Levy's. 

In recent papers, Williams (refs. 4 and 5 ) has outlined a method for 
thick-skinned thin wings with closely spaced spars and ribs (the effec- 
tiveness of which is included in the skin) . In this method, the partial- 
differential equations of plate theory are assumed to be applicable . The 
differential equations are replaced by difference equations in which the 
unknowns are the deflections at a large number of lattice points spread 
over the wing. The method is very much like a relaxation method except 
the "relaxing" is done all at once and for any load distribution by means 
of a matrix inversion. 

Other methods of a quite different nature exist for analyzing low- 
aspect-ratio wings, namely analog methods. These range anywhere from 
testing a model of the wing to building an electrical analog of an ideal- 
ized wing structure as in the method of MacNeal and Benscoter (ref. 6). 

The method to be described in the present paper is an analytical one 
and is essentially an extension to built-up wings of a method developed 
for thin solid wings by Reissner and Stein in reference 7 and "by Stein, 
Anderson, and Hedgepeth in reference 8. The chief difference between the 
theories of Levy, Schuerch, and Williams and the present theory is that 
they simplify the problem by idealizing the structure whereas the approach 
of the present paper is to simplify the problem by idealizing the defor- 
mations . The application of the method is therefore not limited to any 
specific range of skin thicknesses. 

The plate -like nature of the thin low -aspect-ratio wing suggests the 
adoption of some of the assumptions of plate theory. Accordingly, in 
this paper, the displacements of the material points of the wing are 
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assumed to be expressible in terms of tbe vertical displacements of a 
neutral surface. (See fig. 3*) The material points on a normal to the 
undeformed neutral surface are assumed to remain on that normal during 
deformation of the wing. Furthermore , the deflections of the neutral 
surface are assumed to be either linear or parabolic in the chordwise 
direction (both theories are treated in this report), but are otherwise 
unrestricted. An appropriate set of equilibrium equations is obtained 
by using the foregoing specification of displacements in conjunction with 
the principle of minimum potential energy. By approximating the energy 
integral by a finite sum in which all derivatives in the integrand are 
expressed in difference form, the minimization process leads to equilib- 
rium equations in matrix form. 

Specific directions for setting up the required matrices from the 
"raw" data (for example, from drawings of the wing) are given. Two 
load-support conditions are considered fundamental, namely, symmetric 
and antisynmetric . In each of these two cases, a set of influence 
coefficients is derived for a wing that is assumed to be supported in 
a particular way specifically chosen to facilitate routine application 
of the method. (See fig. 4.) A procedure is then given whereby these 
influence coefficients may be modified to obtain those appropriate to 
many other kinds of support including pin- jointed mounting on a flexible 
fuselage . 


SYMBOLS 


pco(y) 


stiffness coefficient for cover plates. 


Dx^dx 


UCq(y) 

Ci(y), Cg(y) functions defining plan form (see fig. 5) 

D local flexural stiffness of wing given by equation ( 5 ) 

E modulus of elasticity of material 

I moment of inertia of spar or rib 

integers 


J; k, 1 

m, n, r, sj 


semispan 

number assigned to station at tip of wing, 1 = Ne 
load at mth reference point on chord at nth station 
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P 

t 

t u> 


¥. 


n,m 


w 


x, y, z 
x g (y) 

y r 


transverse load per unit area 
tMckness of a cover plate 

thicknesses of upper and lower .cover plates , respectively 

deflection at mth reference point on chord at nth station 

lateral deflection of wing (see fig. 5) 
coordinates defined in figure 5 

function defining location of sth spar or stringer (see fig. 5) 
spanwise location of rth rib or stiffener 


z u> z l 


0 


?r 


e 

e x> € y> 


distance from neutral surface to upper and lower cover plates, 
respectively 

spar or stringer stiffness coefficient 
rth rib or stiffener stiffness coefficient 

strain in stringer; distance between equally spaced stations 
7 xy components of normal and shearing strain in cover plates 


£, 9 coefficients of rigid body translation and pitch, respectively 

A angle of sweep of spar or stringer (see fig. 5) 

(i Poisson's ratio 

! n m distance from y-axis to mth reference point on chord at 

nth station (see fig. 7) 

cp k (y) deflection coefficients, (see eq. (l)) 

II energy 
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METHOD OF ANALYSIS 


The structure considered in this paper is a thin wing of low aspect 
ratio built up of cover sheet, spars, ribs, and stiffeners. The carry- 
through structure is considered to be part of the wing. The wing is 
supposed to be internally stiffened so that transverse shear deflections 
can be neglected. It will be assumed that there is a neutral surface 
which does not stretch (or at least that its stretching is negligible). 
Then, since tr ans verse shears are neglected, normals to this neutral sur- 
face will remain no rmal during bending of the wing. Under these condi- 
tions, the displacements of all material points in the wing are expressible 
in terms of the lateral deflections of the neutral surface, assumed to be 
given by the equation (see fig. 3) 


v = cp 0 (y) + *?i(y) + 


( 1 ) 


The potential-energy function which is to be minimized is the dif- 
ference between the strain energy and twice the work of the external 
loads . But before the minimization process is carried out the potential- 
energy function is expressed in discrete form by writing all derivatives 
as differences and all integrations as numerical integrations (by the use 
of the trapezoidal rule or suitable variations thereof). The unknowns in 
the expression then become the values of cp^, cp^, and cp^ at a number 

of stations along the span. After properly accounting for geometric 
boundary conditions, minimization of the potential-energy function with 
respect to each of the unknowns leads to a system of simultaneous algebraic 
equations which can be written in matrix form. The unknowns in this system 
of equations are the qp^ values at the station points; terms on the right- 

hand sides of these equations are derived from the loading. The system is 
solved by a matrix inversion. 


Strain- and Potential -Energy Expressions 

The first step in the procedure is the calculation of the strain 
energy of the wing in terms of the deflection of the neutral surface. 

The energies of the various components of the wing will be considered 
separately . 

Covers . - The strain energy of s fetching of one cover plate (say the 
upper) is given by 
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1 

2 





2 < 1 - 



dx dy 


( 2 ) 


The coordinate system chosen in the plan form of the wing is shown in 
figure 5 together with the definitions of C]_, c 2 , and l » The thick- 
ness of the upper cover plate tu may he a function of x and. y. 
According to the fundamental assumptions , the middle -surface strains in 
the cover are given in terms of the deflection w of the neutral sur- 
face by 


e 


x ~ 



a y 


" " Zu v 


&L 


7 xy = “ 2Zu Sx By 


(3) 


where ■ ^(x^y) is the distance from the neutral surface of the wing 
to the middle surface of the upper cover plate (z u is assumed to be a 
slowly varying function of x and y) . In terms of w, the strain energy 
of stretching of both cover plates is 


H c 




2(1 



s 2 w afw 
Bx 2 By 2 


A 


dx dy 


(4) 


where 


D(x,y) = — l__Lz u 2 + fyzA 
1 — ' ' 


(5) 


in which the subscripts u and l refer to the upper and lower covers, 
respectively. The strain energy of bending of the cover plates can be 
neglected if the thickness of the cover plates is small in comparison 
with the thickness of the wing. This simplification has been assumed 
here. 


When the expression for w given by equation (l) is introduced into 
equation (4), the integration with respect to x may be carried out; the 
result is 
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n c = jao^o") 2 + 2a 1 'Pb" < Pi" + a 2 [(H "f + 2c Po'V'] + 

2a^ L ,, q^" + a^(cp 2 ") 2 + 4aQcp 2 2 + ^(a^cpo" + a^" + a 2 cp 2 ")cp 2 + 

2(1 - p.) jao^i') 2 + ^a^'cpg’ + 4a 2 (cpg ’) j- dy (6) 

where the primes indicate differentiation with respect to y and 


p c 2 

%(y) = / ED^dx 

Jc l 


(7) 


Spars and stringers .- The strain energy of a spar (or stringer) 
running from y 0 to y-]_ at an angle A to the y-axis is 


i r p(yi) nfoi^ao - i r yi 

2 ^p(y 0 ) W 2 / ^ 2 yo 


where p is distance along the spar. 




( 8 ) 


0 = El(y) cos^A 


(9) 


and w(y) = w[x(y),y] where x = x(y) is the equation of the line 
along which the spar lies (see fig. 5)* The moment of inertia I is 
calculated for sections normal to the lengthwise direction of the span 
about an axis lying in the neutral surface of the wing. In terms of 
cpQj <P]_> and cp 2 , the strain-energy expression for a spar (or stringer) 

is 


> y l 


n s = if Ps[<Po” + (VPl)" + ( x s 2c P 2 )”] 2d y 
y 0 


( 10 ) 


where the subscript s is used to identify the spar or stringer. 



8 


NACA TN 5640 


Ribs and chordwise stiffeners . - The strain, energy of a rib (or 
chord-wise stiffener) located at y = y r is 



where I r (x) is the moment of inertia of. ubhe rib about the neutral 
surface. In terms of cpgj the expression for H r is 


where 





Ely dx 


( 12 ) 


( 13 ) 


Loads . - The potential-energy function of the transverse loads of 
Intensity p is 


IXp = -jT l J 2 pw dx dy 


( 14 ) 


A more explicit form for Up when the loads are concentrated is 
given later . 


Finite -Difference "Forms 

The potential-energy function which is to be minimized to obtain 
appropriate equilibrium equations is 


n - n c + ^ n s + ^ -^r + Up (15) 

s r 

An application of the calculus of variations at this point would lead to 
differential equations of equilibrium which could then be put in difference 
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VI 


form.; however. It is easier to proceed to these directly. The minimiza- 
tion process leads directly to a set of simultaneous linear algebraic 
equations for the <p^ values at discrete points along the span when the 

t derivatives occurring in the expression for U are replaced by finite 

differences and the integrations are replaced by finite sums . In this 
paper, the stations along the span of the wing are equally spaced and 
numbered from n = 0 at the center line to n = N at the tip. There 
are also stations corresponding tc n = -1 and n = N + 1. (See fig. 6.) 
The following difference approximations for first and second derivatives 
are vised: 



where n indicates the station and e is the distance between stations . 
Integrals are approximated by finite sums according to the trapezoidal 
rule or, if necessary, by some appropriate modification thereof. Thus, 
the discrete form of the contribution to H c (see eq. (6)) due to a 

typical term involving a second derivative is 


1 

2 



^(■Po") 2 ^ - Uy 


! a O,o('PO,-J 


2 <Po,0 + c Po,i)‘ 


+ 


a 0,l(?0,0 - ^0,1 + *0,2 f + • • • + 

^N-l^,!*-^ “ 2 c ?0,N-1 + c P0,n) 2 + 

2 S OjN( c PO ,N-1 “ 2c P0,N + c ?o,n+i) 2 | 


(it) 


-4 
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or due to a typical term involving a first derivative is 



2(1 - |i)a 0 ^ 1 ') 2 dy = 




+ 


i 


a 0,N- 



(18) 


The discrete form of the other terms entering into the expression 
for H c may he written in a similar way; they are omitted here for 

simplicity. These equations are valid if the a^ values are continuous 
functions of y. However, discontinuities due to cutouts, reinforcement, 
and so forth may occur, in which case some provision for modifying the 
trapezoidal rule to improve its accuracy when applied to discontinuous 
functions should be made. This correction can be made by replacing the 
values of a^ at the station nearest the discontinuity by modified 

a^ values according to a simple rule which will be given later. For a 
typical spar which begins, for example, at station Hq and ends at sta- 
tion n^, the energy may be written in the following way. First let 


„ »• 9n « + „ + x v, 

s,n OjH s jXi 1 j n SjH 2^n 


( 19 ) 


^s,n = x s,n-l ” 2X s,n + X s,n+1 
Then, from equation (10), 

118 = 2€^^ S,ri 0 S}1 ^ 2 + 


+ • • . + 


( 20 ) 


S ,i!r , 2 + ^=6 \|r 2 ] 

K s,n-|_-l r s,n 1 -l 2^s,n]_ Y s,ni J 


( 21 ) 
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which may "be written 


- ^(fs,0+s,0 2 + • • • + h,*s,/) < 22 > 

where the definition of in terms of p g is obvious by a comparison 

of equations (21) and (22). In case the spar does not end exactly at 
station np (or n^), the value of p g at station n^ (or n^) must be 

redefined to take this into account; the exact definition is given later. 
In any case, formula (22) will still hold if f3 SjI1 is properly defined. 

For a rib at y r = (n^ + d^e, the energy H r is 


Hr = *7. 


2 r 


4 - + Vf 2 , v i 


(23) 


where 0 ^ dj. ^ 1. Here linear interpolation has been used to obtain a 
value of <p 2 between stations. 


Boundary Conditions 

The energy expressions in the last section Involve cp^ values at 

stations -1 and N + 1; the disposition of these terms depends on the bound- 
ary conditions. The support system should be chosen so that the boundary 
conditions take a simple form and so that the generality of the result 
is not restricted; that Is, the influence coefficients appropriate to the 
chosen support system must be such that influence coefficients appropriate 
to any other support system can be derived from them (as is not the case 
for cantilever support, for example). For deriving a set of Influence 
coefficients for symmetric loading, the wing is considered to be clamped at 
the origin. (See fig. 4(a).) The conditions on are 


*0 = 


*1 = 


= 'Pi “ 


= 0 


I 


I 


(at y = 0) (24) 
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or in difference form 


^ 0,0 ~ $ 1,0 ~ 0 ^ 0,-1 " ^ 0,1 ^ 1,-1 “ ^ 1,1 ^ 2,-1 “ * 2,1 

(25) 

In the antisymmetric case the wing is considered to he simply supported 
at the points (x=0,y=±e) and constrained (reactionlessly) so that the 
center line does not deflect. (See fig. 4(b).) In difference form 


^ 0,0 ~ * 1,0 “ * 2,0 = * 0 ,- 1 " * 0,1 " 0 


* 1,-1 “ “* 1,1 * 2,-1 _ "* 2,1 

(26) 


At the tip there are no geometrical 
tions on the values of cp^ ^ or cp^. 


contraints and hence no restric- 
to be enforced. The equations 


an 

^*£, 11+1 


= 0 


( 2 ?) 


could be ■used to eliminate the three unknowns cp^ from the energy 

expression, but it is much simpler to retain them as unknowns. 


Matrix Form of Equilibrium Equations 

The final form of the equilibrium equations obtained by minimizing 
the potential -energy function IT subject to the boundary conditions will 
be given in this section. The derivations of the results are mostly a 
matter of routine but lengthy algebra, and they have been omitted from the 
paper except for an example derivation given in appendix B. 

The equations 


an 


= 0 


(28) 
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lead directly to the matrix equation 


A ] [cp] = c 3 W 


or in expanded form 


poo A 01 ^02] K 


*10 *11 *12 *1 = e5 Pi (: 

A 20 *21 *22 9 2 P 2 

The definitions of the matrices [a^-] are given as follows (in which 
primes indicate the transpose of a matrix) : 

M ■ W ([*»] * z [®«])[«o] 1 
W - W(HW’*i:[ 5 j[x.]) 

M ■ M (H H ' * e W M ' ♦ *44) 

W ■ WH H ' + z N H [*.]' * 2(1 - “PI M H ' 

N ■ H H N ' * z H M M’ * *a - >*>. 2 [■>,] ^ ^ ■ * *4^] 

M ■ W W W ' * Z M [s.j [*. 2 ] ' * 8(1 . „ )c2 |^ ■ , 

2,1,2 jw h] * (N w)|+ * 4 >) + « 3 Z [rj 

N-N’ N-W' H-M 1 
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The def ini tions of the matrices entering into the right-hand sides 
of equations (31) and the [qpj and [p] matrices in equation (30) depend 
on the boundary conditions. They are as follows (where the number of 
rows and columns, in that order, are indicated below the equation number): 


For the symmetric case: 

2-2 1 
1-2 1 
1-2 1 


M - H - 


L 


1 -2 1 
1 -2 

j 

1 


(32) 

(N+l)x(N+l) 


-2 1 
2-2 1 



1-2 1 
1-2 1 

1-2 1 
1 -2 

1 


* 


(33) 

(N+2)x(N+l) 


1 -1 



1 -1 


1 -1 


1 


(34) 

(N+l )xN 


u 


0 
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-1 


1 -1 


1 -1 



1 -1 


1 

0 


(35) 

(H+2)xN 


2®*,° 


^k, 1 


a k,2 



a k,N-l 

jpk,!! 


(36) 

(N+l)x(W+l) 


^kjO 


®k,l 


a k,2 



4 


a k,N-l 


ia 

2 


k,N 


0 


( 3T ) 

(N+l)x(N+2) 
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(58) 

NxN 


^ 0,0 




* 0,2 



^^N-l 

2®0,N 


OJ 


(59) 

(N+2)x(N+2) 


* 


Ps,0 




m 

(N+l)x(N+l) 


W- 



M 
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M 


*a,-l + *s,l -2*8,1 *s,l 

*s,2 _2 x s,2 *s,2 

*8,3 - ac s,3 *8,3 


*s,H-l - 23 c b,H-1 *8,11-1 
*s,n - 2x s,ir 
*8,11+1 


M 


- 2 * 8,0 
2 , „ 2 


*s,0 2 



-2*s,l 2 

*s,l 2 


_ 2 
*8,2 

-2*8, 2 2 

*s,2 2 


*8,3 2 

-2x s,3 2 


*s,N-i a -2 *s,H- 1 2 


When the rth rih falls between stations, 




*a,N-l 


*B,N -2*8,^ 




(*L) 

(H+l)x(N+l) 


(te) 

(H+2)x(H+l) 


has the following form: 


0 



r,(i - V 2 ) 

7 r < Jr(l - dr ) 



(43) 

(N+2)x(N+2) 


0 


0 
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Rules for proper positioning of the elements in the |^r r J matrix are given 

in the section entitled "Mechanics of Application." The expanded form of 
the [cpj matrices occurring in equation ( 30 ) is 


* 0,1 

* 0,2 



(44a) 

(N+1)X1 


*0,N+lj 


* 1,1 

* 1,2 



(44b) 

(N+l)xl 


*1,N+1 




9. 


2,0 


9, 


2,1 


*2,N+1 


(44c) 

(N+2)xl 
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Also 



P 0,2 



(45a) 

(N+1)X1 



p iA 

p l,2 



( 45 ^) 

(N+l)xl 





(45c) 

(N+2)X1 


4 


p 2,N 

_ 0 _ 
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where 


Slip 


For the antisymmetric case: 


( 46 ) 


M - 


1-2 1 
1-2 1 


• • • 


1-2 1 
1 -2 

1 


( 47 ) 

NXN 


D 1 


n 


,] - 


-2 1 
1-2 1 


i i « 


1 -2 
1 


1 

-2 

1 


( 48 ) 

(N+1)XN 


N ■ W - 


i -i 
i -l 


i -i 
i 


(^9) 

(N+1)XN 


0 
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a k,l 

a k ,2 



a k,N-l 



( 50 ) 

NXN 





a k,2 


®k,N-l 

IK,H 0 


( 51 ) 

Nx(N+l) 



same as symmetric case 


® 0,1 



a 0,2 


^N-l 

i®0,N 


( 52 ) 

(N+l)x(N+l) 


0 
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Ps,l 



■ 2x s,l x s,1 


x s ,2 -2x B)2 

x s ,2 

*a,3 





x s,N-l _2 x s,N-1 x s,N-X 


-2x s,K 

x s,N+l 


~ 2X V,\ Z 

X B,1 2 


X 8,2 2 

_2x S,2 2 



V 2 

X B,3 

-2Xs,l S x 8,3 2 


Xs 2 - 


Y 2 Oy 2 Y 2 

^N-l _2 x s,N-1 *8 ,N-1 

x S ,N 2 -2 xs,N 2 

_ 2 

x s,N+l 


( 53 ) 

NXN 


( 54 ) 

(K+l)xN 


(55) 

(N+l)XEf 


The matrix JJT r j has a form similar to that of the symmetric case except 
that it is of order (N+l)x(N+l). The proper placement of the elements 
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in |^r r j will be given in the section entitled "Mechanics of Applica- 
tion. " The expanded form of the [p^ matrices occurring in 
equation ( 30 ) for the antisymmetric case is 


^ 0,2 

90,3 



( 56 a) 

NXl 


9q,N+1 


91.1 

91.2 



( 56 h) 

(N+l)xl 


9l,N+l 


9 2 .1 

9 2 .2 



( 56 c) 

(N+l)xl 


*L 


92,N+1 
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p 0,2 

p o,3 



(57a) 

Nxl 


P 0,N 

0 


p i,i 

p l,2 



(57b) 

(N+1)X1 


P 1,N 
. 0 . 


p 2,l 


p 2 ; 2 



(57c) 

(N+1)X1 


P 2,N 


0 
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Inversion of J~a] Matrix 


The next step in the derivation is the solution of the system of 
equation (30). This may be done either by inverting [A.J at once or 
alternatively by first eliminating jq>oj (as in ref. 8) and then 

inverting a matrix of roughly two-thirds the order of £a] . The most 

economical choice depends on the number of stations and the computing 
facilities available. 


The result of eliminating 



is 


*12. ' A io A oo 1a 01 ! 

I 

A 1 2 ~ A 10 A 00 A 02 



c 3 

P 1 “ A lcAx) ^0 

r“l 

H-p 

1 

0 

-g 

1 

« — 1 

A 22 " ^20^00 ^2 



= 6^ 

P2 " A 20 A 00 _1 PO 


(58) 

or say 


B n B i 2 


<Pl 

^21 B 22 


qp 2 



PO 

Pi 


P 2 


(59) 


where fxl is the identity matrix. There would hardly be any advantage 

' — * i- -i ~\ 


in eliminating 


N 


if setting up the 


M 


matrices for computation 


were as complicated a task as it looks . However, the terms which involve 


products of 


A q 


matrices can be considerably simplified and lead to 


the following results: 
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H ■ W - (H H + ? H hj) (H + ? Hf (H H + ? H W) 

H = W • (H H hlhj (H + 5 hi) (jhlH h 2 J M + ^N’) 

H ” W ■ (h] H + ^ M M + ^ 2 N ) (W + 5hl) (W H hi W + ') 

[®ai] = [^m] 


Note that there are essentially only three different matrices within 
parentheses and that one is the inverse of a diagonal matrix. Also 


(60) 


- (K 




n u y (hoi + n [pb] 

S L J L jy \L J s L J 

II ^s 2 fs_ + ^ W) ([“o] 


1 r VI 


-1 r i-i 

d 0 


(61) 


where in the symmetric case 


i 1 


M 


-1 


2 2 1 
2 2 


12 5^ 
12 5 


1 2 
1 


N + 1 
2 

N 

N - 1 

5 

2 


(62) 


1 



NACA ON 3640 


27 


and in the antisymmetric case 


12 3 4 



12 3 


N - 1 



12 3 

1 2 

1 


(63) 


In order to complete the solution for 
matrix which must be inverted is 


[cp] in terms of PpJ , the 



B. 


'll 


B, 


21 


B 22 


( 64 ) 


The solution may then he written 


_ _ 


f 

— » — 1 


- 


1 

r "i 

*0 


1 

(C'B _1 )C 1 C'B -1 


^o' 1 ! 0 

"T "to" 
1 
1 

0 ! 0 

0 

\ 

Po 

q>i 

= e5 


(c’B -1 )'! B " 1 

+ 

0 


Pi 

cp 2 


\ 

\ ' | 

1 


0 

- 

/ 

F 2 . 


or 



where 



( 63 ) 


(66) 


(67) 
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and 


>] _1 ■ ([ D o ] _ 1 )'([ a o ] ps ])" 1 ^]' 1 <68) 


In certain cases, the preceding operations will have to be modified 
because the matrices which are to be inverted turn out to be singular. In 
particular, this may happen when the wing has a pointed tip. The reason 
the matrices are singular is that at some stations the three unknown sets 
of cpk values are linearly dependent; in other words, the deflection is 

over determined at that station. The details of the remedy for this situa- 
tion are given in the section "Mechanics of Application." 


After the [a]"^ matrix is obtained, it is modified slightly. 
According to the derivation there are rows and columns of £aJ“^ corre- 
sponding to cp^ and p^ These three rows and columns are to be 

deleted. Also, in the symmetric case, it will prove convenient to insert 
two rows of zeros corresponding to cpQ^Q and cp-j_ Q and ^ w0 colu™ 18 0 f 

zeros corresponding to pq q and pj_ g. In the antisymmetric case, 
insert a row of zeros corresponding to cp^ ^ and a column of zeros corre- 
sponding to Pq meaning of jjpj and JjpJ is altered but the 

notation will be kept the same In what follows. Also, the factor 


(see eq. (5°)) is included in the modified [a] - 1 matrix. The result is 
a symmetric matrix which will be called [jg] . Equation (66) becomes now 



(69) 


Equation (69) gives generalized deflections In terms of generalized loads 
and essentially expresses the solution to the load-deflection problem. 

For some applications (finding modes and frequencies, for instance), it is 
convenient to proceed from this point keeping the cp^ as generalized 

deflections . Only the problem of finding influence coefficients will be 
worked out In detail. 


Influence Coefficients 

A set of reference points on the plan form of the wing must be chosen 
with respect to which the matrix of influence coefficients will be given. 
The reference points on the plan form of the wing are located in terms of 
the set of numbers £ n m which give the distance from the y-axis to the 
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mth reference point on the chord at the nth. spanwise station. (See 
fig . 7 . ) A quantity referring to a particular reference point will "be 
labeled with the subscripts n,m. There may be any number (cl^J of ref- 
erence points along the chord at station n but the deflections at only 
three of the reference points are independent in the parabolic theory 
(two in linear theory). 

The relation of JjpJ to the concentrated loads P n ^ m at the refer- 
ence points can be obtained through the definition of p^ n in terms of 
Up; that is. 


P 


k,n 



(70) 


In terms of _ and W_ _ (the deflection at reference point (n,m)), 

ll ^ ILL IJ. y ill 

the expression for lip is 

Up = w n,m?n,m (T3-) 

n,m 
or 

^p = ^ (^0,n' + ^n^m^l^n + ^nF^nVn,]!! (f2) 

n.m 


Then from equation ( 70 ), 


p k,n * H 
m 



(73) 


which may be written in matrix form as 



( 74 ) 


or in expanded form, for the symmetric case, 
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In the antisymmetric case, [h] is given by 
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Fran the definition of w in terms of cp^ (eq. (l)), the following 
matrix relation may be deduced: 


H-H'W (77) 

where JVj is a column matrix whose elements W n m are arranged in the 
same order as the elements P n m in [p] . (See eq.. (75).) From equa- 
tions ( 69 ), (74), and ( 77 ), it is evident that 


where 



(78) 



(79) 


The matrix [g] gives the desired set of influence coefficients. 

As defined in this paper, the influence coefficients given by equa- 
tion ( 79 ) are the deflections at symmetrically placed reference points 
(with respect to the center line) due to symmetrically placed unit loads. 
There is one instance in which the definition may be a little ambiguous 
and that is when the load is on the center line in the symmetric case. 

In this case the deflections at all reference points are to be interpreted 
as due to a double load on the center line since now the two unit loads, 
symmetrically placed, have moved into coincidence. The obvious alternate 
definition would not result in a symmetric [g! matrix. 


Equations of the Linear Theory 


The equations of the linear theory may be obtained from those of 
the parabolic theory by setting ^2 n = 0 and - omitting the equations 
^ | — 1‘ 

= 0 at the outset. The steps in the subsequent derivation are 


K 2 ,n 
identical. 


In this case the matrix which has to be inverted is either 


Ado 



Adi 

A n 

_ 


(80a) 



32 


NACA TN 5640 


or 

M - [ E n] 

In the latter alternative (that is, when [b] is inverted^, 
be obtained from the equation 


(80b) 

[a] -1 can 


where now 


C’B" 3 ^ 

1 " 

i C’B " 1 
{ 


1 1 

1 

1 

0 

! 

H 1 

1 1 

Ol 

lW 

(c’B" 1 ) 

1 B- 1 

1 

+ 

0 ! 0 

1 


[C] . [cj 


(81) 


(82) 


The same modification of [a] - 1 as is described in the paragraph preceding 
equation '(69) is again made to obtain jj^|. As before 


[°] - WWW (85) 

where j is the same as in parabolic theory except the last , N rows 
are omitted. 


Adaptation to Other Loading or Support Conditions 

A knowledge of the load-deflection characteristics of the wing apart 
from the fuselage is often insufficient for accurate aeroelastic analysis 
of the airplane. However, the effect of attaching the wing (which 
includes the carry-through structure) to the fuselage at a number of sim- 
ple supports (no moments transmitted) can be accounted for by a fairly 
simple modification of the influence coefficients obtained by using the 
standardized supports . The details of the analysis are given in 
appendix A. 

In the analysis of the present paper, the carry-through structure 
has been assumed to be part of the wing. The assumption is valid if 
the carry -through structure is not very different from the wing struc- 
ture; for example, this would be the case If one cover sheet and some 
stringers were cut out but the spars were continued on through. If 



NACA TN 3640 


53 


the carry-through structure is appreciably different from the wing struc- 
ture, then a modified, analysis is called for. The analysis of the present 
paper could probably be adapted to the case in which each half -wing is 
attached to a flexible structure characterized by a set of influence coef- 
ficients which involve rotations and couples as well as deflections and 
forces at the points of attachment of the half -wings . 

For some applications it is necessary to know the deflections due to 
a unit couple applied at some point along the trailing edge (as caused by 
an aileron, for example). If the trailing edge is parallel to the y-axis 
and the point of application of the couple falls on a station, then the 
solution to the problem is simple. A reference point is introduced off 
the trailing edge at the reflection of some reference point on the chord. 
Equal and opposite loads at these two reference points produce the same 
effect as a couple applied at the trailing edge within the framework of 
the parabolic (or linear) theory. The deflection due to a load at the 
reference point off the wing is found in exactly the same way as though 
it were on the wing. In case the trailing edge is swept or the couple is 
applied between stations, the correct procedure for finding the deflections 
due to the couple can be found by recourse to the energy method (find Up)* 


DISCUSSION 

Approximations Involved in the Theory 


The approximations involved in the present theory are, for the most 
part, consequences of restrictions placed on the displacements at the 
outset. Seme of these approximations, such as the neglect of transverse 
shear deflections, are quite common and need little discussion but seme 
others, more or less peculiar to the present theory, require some words 
of justification and some warning as to the limitations they place on the 
theory. 

The restriction to parabolic (or linear) chordwise deflections has 
already been used subcessfully in the case of a solid plate (ref. 8). In 
the present case, the ribs near the root are not expected to bend appre- 
ciably because of the restraint offered by the fuselage. For a wing of 
low aspect ratio, nearly the whole wing is within the region of influence 
of the root restraint. Also, since the ribs have zero (or, at most, small) 
bending moment at either end, the chordwise bending moments are much less 
than the spanwise bending moments so that, if the stiffnesses in the two 
directions are of the same order, the chordwise curvature should be much 
less than the spanwise curvature. Thus, provided the ribs are not too 
light, it seems justifiable to assume that the chordwise deflections may 
be well approximated by a parabola, or even by a straight line. 
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If the wing is perfectly symmetrical, top and. bottom, so that the 
middle plane may be taken as the neutral surface, then the stretching of ‘ 
the "neutral" surface is certainly extremely small (zero within the 
framework of small-deflection theory) when the wing is bent under lat- 
eral loads. However, if there are cutouts on the top or bottom of the 
wing or various other asymmetries, then the situation is not nearly so 
clean-cut. The difficulty is that there does not necessarily exist a 
surface which remains unstretched (even in small -deflection theory) when 
the wing is subjected to a lateral load. As an exaggerated example, 
suppose there are many spanwise stiffeners on the top of the wing and 
many chordwise stiffeners on the bottom. In such a case, it would seem 
reasonable to specify the u and v displacements with respect to two 
surfaces, one on which u is zero and the other on which v is zero, 
rather than with respect to a single neutral surface on which both u 
and v are zero. That such a case would arise in practice is unlikely, 
but at least it shows the need for making some reasonable assumptions. 

The only way really to avoid the necessity of making engineering judg- 
ments with respect to the choice of a neutral surface would be to go to 
some more exact theory in which more complicated expressions for speci- 
fying displacements are assumed. Barring that, the locus of the principal 
axis of inertia of the chordwise cross sections seems to be a reasonable 
choice. The resulting influence coefficients probably will be in error in 
the neighborhood of a cutout or other discontinuity, but the overall 
results should be accurate. 

The method of analysis of the present paper has been motivated by the 
resemblance between the low-aspect-ratio wing and the plate inasmuch as 
they are both thin and flat. However, since the carry-through structure 
may not even remotely resemble a plate, the validity of the method must be 
examined on that point. The criterion is whether the assumed form of the 
displacements can accurately describe the true shape of the deformed struc- 
ture. If, for -example, the carry-through^ structure is composed of three 
beams which are the continuation of three wing spars, then the deflections 
of these three beamB can be described exactly in terms of the three func- 
tions cp 0 , cp^, and cp^ of the parabolic theory. In this case, there is 

no question that the analysis is adequate as far as the carry-through 
member is concerned.' However, the chordwise location of the continuous 
spars could conceivably be such that the deflected shapes of the Inboard 
chords qf the wing could not be approximated very well by a parabola even 
when the loading on the wing is reasonably well distributed. In addition, 
the discontinuities in the structure cause discontinuities in the higher 
derivatives of cp^ and this may cause an appreciable error in the differ- 
ence equivalents of the lower derivatives used in the analysis. 
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Refinements Over Beam Theory 


Some of the ways in which the present theory falls short of exact- 
ness were pointed out in the last section. There' are, in favor of the 
present theory, a number of refinements over beam theory which are either 
mani fest or implicit in the equations . The additional generality of the 
more refined theory allows for restraint against warping of the root 
cross section and for restraint against anticlastic curvature at the 
root. In a long beam these effects are felt only in the region near the 
root, but in a low-aspect-ratio wing where the root chord may be about 
as large as the semispan these effects are appreciable over the whole 
wing. A glance at the form of the equations shows that bending ^cpQ j , 

torsion ^qp^), and chordwise bending ^p 2 ) are inextricably coupled 


together; however, in simple beam theory, the equations are uncoupled 
by virtue of the assumption of an elastic axis . 


Comparison With Other Theories 

In either Levy's or Schuerch's theory the direct-stress -carrying 
capacity and the shear-carrying capacity of the sheet are separated. 

The direct-stress -carrying capacity of the sheet 1 b lumped in with the 
spars and ribs. This is a common assumption in the' treatment of semi- 
monocoque structures with thin skin. However, when the cover plates get 
thicker, and thus form a proportionately greater part of the total wing 
material, a more refined treatment of them is appropriate. Probably the 
most important effect neglected in the thin-sheet approximation is the 
coupling between spanwise and chordwise stretching of the sheet - in 
other words, the Poisson's ratio effect which produces anticlastic 
curvature. In the present theory, a nearly exact expression for the 
strain energy of the cover sheets is used. The accuracy of the result 
depends on how well the assumed displacements approximate the true ones, 
but does not depend directly on the thinness of the cover sheets. So, 
at least in*the parabolic theory, the Poisson's ratio effect is taken 
into account. Williams' method is the antithesis of Levy's in the sense 
that primary consideration is given to the covers rather than the spars 
and ribs . However, Williams 1 method is not well adapted to the analysis 
of a wing with a few heavy spars and ribs or with thin skin. At the 
present time the ranges of applicability of the several theories have 
not been well defined owing partly to a lack of experimental evidence and 
partly to the newness of the problem. 


Some Possible Improvements and Extensions of the Theory 

The effect of transverse shear deflections in the ribs and spars is 
possibly the most serious omission in the present theory as it now stands.- 
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For some configurations the effect of shear deflection in the ribs 
would be more imp ortant than bending of the ribs. Experience may indi- 
cate that for practical designs appreciable error is made by neglecting 
shear deflections, but the extension of the theory to take this into 
account has been left for possible future work. (See, however, ref. 5*) 

The present method could be generalized to allow arbitrary chord- 
wise bending. The generalization would involve replacing the double 
integral for the strain energy of the cover sheet by a double sum. The 
unknowns in the equations would be the values of w at many lattice 
points on the wing rather than the values of at a number of span- 

wise stations. The resulting theory would be similar to Williams' theory 
except that the ribs and spars would not be spread out to act with the 
sheet . 

Some more precise treatment of chordwise bending would be necessary 
for a complete stress analysis of delta wings. The primary concern of this 
paper has been deflection analysis and, although the stresses can be esti- 
mated from the deflections, a double numerical differentiation is involved; 
consequently, the resulting stresses would be less accurate than the 
deflections. The parabolic or linear theory would give either constant 
or zero curvatures of the ribs; thus, the theory would seem to be of 
questionable value for determining stresses in these members. However, 
in reference 8 a comparison of experimental results with the predictions 
of the parabolic theory applied to a 8 olid delta plate shows good agree- 
ment for the spanwise bending stresses. 


MECHANICS OF APPLICATION 


More theory and mathematical details have been given in the section 
entitled "Method of Analysis" than are necessary in actually setting up 
and solving a given problem. On the other hand, some troublesome points 
which may arise in the application of the method have been omitted so far. 
This section in intended to serve as a guide in applying the method and 
also to cover these troublesome details. 


Steps in Setting up Matrices of Structural Properties 

The following steps are used in setting up matrices of structural 
properties: 

(l) Choose a coordinate system as in figure 5 with the origin on 
the center line. 
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(2) Choose a number of equally spaced stations along the hal f -span 
as in figure 6. There may be some "natural" number of stations dictated 
by the geometry of the structure but, in general, reasonable accuracy 
can be obtained with about eight stations . For convenience, the points 
of attachment of the wing to the fuselage should be made to fall on sta- 
tions, if possible. 


(3) Choose a neutral surface. The choice recommended in this paper 
is the locus of the principal centroidal axis of the chordwise sections. 
However, any discontinuities such as those due to cutouts and reinforce- 
ments should be smoothed out. Of course, in the process of finding the 
neutral surface the cross-section areas of spars and stringers, skin thick- 
nesses, and so forth will be tabulated at the stations and half -stations . 
The effective cross-section areas of swept spars and stringers should be 
taken to be the actual cross-section areas (normal to their own axis) 

times a reduction factor of cos^ A a , where Ag is the angle of sweep of 
the spar or stringer. 


(4) Calculate tables of the coefficients a^^ 
the formulas 


and 


a k,n+i 


from 


a k,n 



DqX^cIx 


(k = 0, 1, 2, 3, *0 
(n =* 0, 1, 2, . . . N) 


(84) 


*k 


1 = r 2 1^,1 x k dx , 

> n+ a J „ p (n - 0, 1, 2, . . . N-l) 


(83) 


where (see eqs. ( 5 ) and ( 7 )) 

D = + tlZl 2 ) 


( 86 ) 


These integrations can be done numerically if necessary. In case there 
eire discontinuities in the cover stiffness D, seme adjustment of the 
values of a^. n at the station nearest the discontinuity is necessary. 

Let the jump in be (positive if a^ increases across the dis- 

continuity for increasing x) and let ed be the absolute value of the 
distance between the discontinuity and the nearest station m; then, 

(a) If the discontinuity in a^ is a distance ed to the 
left of station m, replace a^ m by a^ m ~ (h ~ 
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(t>) If the discontinuity in is a distance €d to the 

right of station m, replace a k ^ m by a^ m + ^ - 

(c) If the discontinuity is at station m, use the mean value. 

A similar rule holds for adjusting the values of the coefficients 
a k 1. It is probably sufficiently accurate to obtain the values of 

j p 

the jumps ^ graphically. 

(5) Calculate the coefficients j3 g n for each spar and stringer at 
each station from the formula 

Ps,n " ( E1 s)n cob5 *s < S 7> 

where I is computed for a cross section perpendicular to the spar about 
an axis lying in the neutral surface and "A s is the angle of sweep of 

the sth spar (stringer); see figure 5* 

(6) Tabulate the coefficients j3 s ^ n * The rules for finding the 
values of (3 in terms of 0 are 

(a) If the left end of a spar (stringer) falls within ie of 

station m, let edg be the signed distance from station m to the 
end of the spar; then, 

3 = (b “ clJ) 3 (88) 

p s,m \2 °j H s,m 

(b ) If the right end of a spar falls within ie of station m’, 

let edg' be the signed distance from station m' to the end of the 
spar; then, 

3 s ,m' = (|+ (®9) 

(c) For m < n < m' , 

3 =3 

p s,n p s,n 


( 90 ) 
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(d) If the spar crosses the center line, 

K, 0 - ¥s,o od 

(e) Elsewhere, 

h,n - 0 (92) 

Sometimes p g will have to be extrapolated to obtain a value at a sta- 
tion point. For a light stringer, little error is committed by assuming 
the stringer to terminate at the nearest half -station on the outboard 
end; in this case, p is the same as p at any station (except at the 
center line) to which the stringer extends. 

(7) Let x s ^ n be the distance from the y-axis to the sth spar 
(stringer) at station n. (See fig. 5.) Tabulate 

(a) x s for the unswept spars and stringers 

(b) x s n for the swept spars (stringers) including values for 
x s,-l ** Ps,0 t 0 and x^ N+1 if p B , N ± 0. In any case, the 
values for x g need only be tabulated for one station beyond the 
ends of the range for which p g ^ 0. Values of x g at stations 
off the span of the half -wing are to be obtained by extrapolation 
( x s,-l = 2*8, 0 " x s, l) * 

(8) Integrate the stiffnesses of the ribs and chordwise stiffeners 
across the chord to obtain a set of values of 7 r according to equa- 
tion (13). For a rib on the center line, take 7 r to be one-half the 

value given by equation (13) • For those ribs or stiffeners that fall 
between stations, the numbers dj, (where ed r is the distance from the 
first station to the left of the rib) should be recorded. 

All the necessary information has now been extracted from the design 
of the wing and arranged in tables of numbers. 

a^j , and pij Beatrices from the values of 

ajj.*J matrices from the values of a^ at 

the half -stations according to equations (36) to (39) (symmetric) or 
(50) to (52) (antisymmetric). 


(9) Set tip the j^a^J , 
a-jj, at the stations and the 



40 


NACA TN 3640 


(10) Set up the [4 [4 and [4 matrices for the swept 

spars and stringers . For spars or stringers not running the full length 
of the span, some of the columns in j3 s will be zero, in which case the 

elements in the corresponding columns of £xjj and can be set 

equal to zero without affecting any subsequent result. For unswept spars 
and stringers, x s = Constant; thus, jjx s j = x s j^l] [x s ^} = Xg^jjjrjJ. 

There is no need for setting up the |jx s j] and [x g ^J matrices in this 


case. A considerable saving in computational labor is afforded by 
including the effects of the unswept spars and stringers in with the 
skin (no approximations involved) in the following way: Set up the 


matrices 


t 


0sx fl 


Hi 4 [rijs4 [Hi ^4 [si^s 5 ]' ^ 




which will be diagonal matrices of the same form as 


Add these matrices directly to [jSqJ, jja^jj , [jag], jja^J , and [jaijj 
respectively ^but not to jja^j , [a^p] , or 


[4 


(11) Set up the matrix [Z r r ] from the table of values of y T 
and dj. where each rib or stiffener contributes terms to the elements in 


the matrix 


[ErJ 


as follows: 


If the rth rib (stiffener) is between stations m and nn-1, a dis- 
tance edj. from station m, then the rth rib contributes the block of 

terms 


y r (l - 4 ^4 “ d r) 

Y T &r(± - d^ 7 r <3r 2 


where the upper left-hand element is in the (nH-l)st row and column in 
the symmetric case and in the mth row and column in the antisymmetric 
case. For ribs at a station, (ij, = 0. There is no contribution from a 

rib on the center line (station 0) in the antisymmetric case, and in 
this case those ribs between the center line and station 1 contribute 


only to the element in the first row and column of 



In either 
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the symmetric or antisymmetric case 


- M 


is a square matrix which, is 


either diagonal or has numbers on the principal diagonal and the two diago- 
nals next to it. In the symmetric case, is an (N+2)x(N+2) matrix 

and in the antisymnetric case it is an (N+l)x(N+l) matrix. In either case, 
the last row and the last column are filled with zeros . 


Calculation of Matrix 

There are two alternative procedures for obtaining the inverse of 
J^aJ . The first alternative involves setting up a matrix of the order of 

about 3NX3N and inverting it. The second alternative involves setting 
up a matrix of the order of about 2NX2N and inverting it and also setting 
up three auxiliary matrices and performing several matrix multiplications 
and additions subsequent to the inversion operation. The most economical 
choice between the two alternatives depends on the computing facilities 
available. In general, the first alternative involves the least nonauto- 
matic computing and should be preferred when good high-speed computing 
equipment is available (provided N is not too large). The second alter- 
native is also suitable for automatic computing, but involves the setting 
up of some additional matrices before the work is ready for the machine. 


First alternative .- Using the results obtained thus far, form the 
matrices |~A_yJ according to equations (31)* Note that the summations 

should be made only over swept spars and stringers if the effects of the 
unswept spars and stringers have been included in the [akj niatrices 
as stated in step (lo) for setting up matrices of structural properties. 
Combine these matrices into the matrix |V| according to equation (30 ). 
The next step is the inversion of the matrix |aJj however, in certain 

cases this matrix will be singular. The proper way in which to modify 
[a] before inversion in certain of these cases is listed as follows: 


(l) In case only two spars (which do not taper to zero) extend to 


a pointed tip, strike out the last row and column of [aJ ^ 
setting <P 2jlJ+1 = oj. 


corresponds to 
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(2) In case only one spar (which does not taper to zero) extends to 
pointed tip, strike out the last row and column in |~a]J and also the row 
and column in [aJ containing the last row and column of [ A ll] ^corre- 
sponds to setting cp^ N+1 = <P 2 ,N +1 = °)* 


(5) In case all spars taper to zero at a pointed tip, strike out 
the rows and columns in £aJ which contain the last rows and columns in 


Aqq > {_ A H > £ A 22 J ^corresponds to_setting cp^ 


<P. 


2 ,W+1 


-o). 


N+l " ^N+l 


(4) If there is a carry-through bay with only two spars and no 
cover sheet or ribs , strike out the rows and columns in ^A ] which con- 
tain the first few rows and columns in • ® ie num 4>er of rows and 

columns deleted is the number of stations within the carry-through bay 
(counting the station at the center line in the symmetric case but not 
counting it in the antisymmetric case ) . This corresponds to setting 
cp 2 0 = ^2 1 = ' ’ * = ^2 m = ^ w ^ ere m is iho last station within the 

carry -through bay. This cane may occur in combination with the first 
three cases. 


Second alternative .- Form the matrices according to equa- 

tions (60) and the [cpj and [c 2 ] matrices according to equations ( 6 l) 
Combine these matrices to form [^B ] given by equation (64) and [^C J 
given by equation ( 67 ). The next step is to obtain Here again, 

however, the matrix may be singular. The appropriate modifications to 
[b] if this is the case are as follows (where the cases are numbered as 
in the first alternative ) : 

(1) Strike out the last row and column in [bJ and the last row 
of [c]. 

(2) Strike out the last row and column in [b 3 and the row and 
column in [jbJ containing the last row and column of jj^llj • Also 

strike out the last row of and j^C 2 ^j. 
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( 3 ) Strike out the last row and column of all matrices 4 appearing 
in equations (60), (6l), and (68) including ^aqj + 

as 'usual. 


M 


and proceed 


( 4 ) Strike out the rows and columns in [b] which contain the first 
few rows and columns in J* Strike out the first few rows in j^CgJ . 

The number of rows or columns to delete is the same as in case ( 4 ) in the 
previous section. 


After either the jjjJ matrix or the modified QbJ matrix has been 

inverted, form the [^a] matrix in equation (66) using the definition 
given by equation (65). 

The effect on the inverse of [~aJ of a modification in the struc- 
ture. - If the [a] matrix has been computed and some change in the 
structure is made that affects only a few of the elements in [~aJ, then 

the work of obtaining the new inverse matrix can be greatly reduced by 
making use of a method given in reference 9 - The new inverse is obtained 

exactly whether or not the changes in the elements of [aT] are small. 

This method would be particularly useful for finding the effect of modi- 
fying one or two ribs, or a spar over a small portion of the span. 


Calculation of the £gj Matrix 
According to equation (66), which is 


H - 


( 93 ) 


there are rows and columns 
Pk N+l* Tf these rows and 


in [a]" 1 corresponding to cp^ ^ ^ 
columns are not already missing in 


and 
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by virtue of their deletion under case (l), (2), (3); or (4) of the 
last section, they must now be eliminated. One way to locate the rows 

in question is to write down the column matrix jjpj. There are as many 
rows in £cpj as there are in j~Aj . Thus, if no rows and columns in [a 
have been deleted, then the rows to be deleted in £a] are in the 
position marked with an asterisk. The corresponding columns in [a ] -1 
are also deleted and the result is a symmetric matrix. 


Symmetrical Antisymmetrical 

case case 


Xi' 

• 


*0,2 

• 

• 

*0,N 


• 

• 

*0,N 

<P 0,N+1 

* 

*0,N+1 

*1,1 

• 


*1,1 

• 

• 

• 

*1,N 


• 

*1,N 

*1,N+1 

* 

*1,N+1 

*2,0 

• 


*2,1 

• 

• 

*2,N 


*2,N 

*2,N+1 

- — 

* 

*2,N+1 

- 
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In addition to the foregoing modifications, two rows and columns of 
zeros should he inserted in the matrix in the positions marked 

with an asterisk in the symmetric case. In the antisymmetric case, 

^Aj - ^" is to he bordered on the top and left hy a row and column of zeros 
as indicated. 


Symmetrical Antisymmetries! 

case case 


0 

* 

0 

*0,1 


*0,2 

• 


• 

• 


• 

• 


. 

*0,N 


*0,N 

0 

* 

*1,1 

*1,1 

• 


• 

• 

• 

• 


*i’n 

*1,N 


*2,1 

*2,0 

• 


* 

• 

• 

*2,N 


• 

CD 

y 2,N 


— 


Finally the modified [a]" 1 matrix is multiplied hy and the 

result called [gj as in equation (69). 


Influence Coefficient Matrix 

Choose a set of reference points on the plan form of the wing for 
which influence coefficients are desired (see fig. 7)* Include among 
them the points of attachment of the wing to the fuselage. These ref- 
erence points must he at the spanwise stations, hut there may he any 
number of them along a chord at a given station. Construct the matrix 

[nl (given hy eq.. ( 75 ) in the symmetric case or hy eq. (76) in the 
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antisymmetric case) from the values of l n m which are the distances 
from the y-axis to the mth reference point on the chord at station n. 
Calculate the influence-coefficient matrix [gJ from equation (79) • 

The influence-coefficient matrix for the wing supported at an arbitrary 
number of flexible supports may be found by applying the method given 
in appendix A. 


CONCLUDING REMARKS 


A method for obtaining influence coefficients for thin low-aspect- 
ratio wings has been presented. The application of the method has been 
organized as much as possible into a routine procedure. 

The technique for arriving at difference equations in matrix form 
directly (rather than via differential equations) has afforded consider- 
able economy of thought particularly in the treatment of boundary condi- 
tions at the wing tip and in the handling of discontinuities. 

The development (after eq. (69)) has been limited to finding influ- 
ence coefficients; however, all problems are not necessarily most con- 
veniently handled by means of influence coefficients. For instance, in 
the problem of finding natural modes and frequencies, it is convenient 
to introduce generalized inertia loads on the right-hand side of equa- 
tion (69) in place of the generalized loads matrix QpQ. The proper 

matrix form for the inertia terms can be deduced from the expression for 
the kinetic energy in discrete form. The frequencies and (generalized) 
modes can be found directly from the matrix equations without any need 
for finding influence coefficients. 

Theoretical influence coefficients obtained by the present method 
have not yet been checked experimentally. Perhaps, for practical pur- 
poses, a close check on influence coefficients is too severe a criterion 
for the usefulness of the theory. That is, there may be greater dis- 
crepancies between predicted and experimental deflections due to a con- 
centrated load than there would be if the load were more uniformly 
distributed. In case of serious discrepancies between predicted and 
experimental results, neglect of the effect of transverse shear deflec- 
tion in the theory should probably be the first assumption to question. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 51 ? 1956 . 
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APPENDIX A 

INFLUENCE COEFFICIENTS FOR THE WING ATTACHED TO THE FUSELAGE 


This appendix will be concerned with the problem of deter min ing a 
set of influence coefficients for the wing which t ake s into- account the 
fact that the wing is attached to a flexible fuselage. These influence 
coefficients will be determined by appropriately modifying the set of 

influence coefficients [gJ appearing in equation (79)* For present 

purposes, the set of influence coefficients [jlj must include those 

coefficients which refer to the points of attachment with the fuselage. 

The only information about the fuselage necessary for the solution of the 
problem is a set of influence coefficients for those points of the fuse- 
lage where the wing is attached. The coefficients for the fuselage need 
not necessarily be obtained from as refined an analysis of the fuselage 
as has been carried out for the wing. For example, they might be obtained 
by approximating the fuselage by a simple beam in the symmetric case or by 
a torsion box in the antisymmetric case, or they may even be taken to be 
zero corresponding to the crude assumption of a rigid fuselage. The influ- 
ence coefficients for the fuselage are supposed to be obtained by assuming 
the fuselage to be supported in such a way as to prevent rigid-body 
motions. However, the reactions at these supports must vanish when the 
fuselage is subjected to any self-equiliorated set of loads . 


Mathematical Derivation in the Symmetric Case 

Let the influence coefficients for the fuselage which refer only 
to the points of attachment of the wing be arranged in the mxm matrix, 

Hs (s for symmetric deformations), where m is the number of attach- 
ment points. If necessary, rearrange the matrix [gJq so that the 

influence coefficients which refer only to the points of attachment appear 
in the upper left-hand corner in the same order as the corresponding coef- 
ficients in [V]g. 


Let the reactions at the m supports be written as the column 


matrix 
column of 


jjR] and let be [rJ augmented by zeros to make up a 


M elements where M is the total number of reference points. 


Let the external loads be written as the column matrix 


Let 


M 


H 


with M rows . 


be a column matrix of m rows whose elements are the deflections 
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at the supports and let W he the column matrix of deflections at all 

M reference points on the wing where a given reference point is at the 
point (£,tj). Furthermore, define the following matrices: 

matrix of the first m rows of _[gJ g 



W 

[i] 

[r] 

to 


column of 
column of 
column of 
column of 


m ones 
M ones 

i's at the first m reference points (supports) 
| ' s at all M reference points 


For symmetric loading the deflection of the wing is given by 

M = [g] s ([p] + Hb) + s[i] + e[e] (A1) 


The last two terms on the right axe the symmetric rigid-body displace- 
ments, namely, vertical translation and pitch. The deflection at the 
supports only is 


[w] = [g] s [p] + [g] s [r] + s[t] + e[r] = -[j] s [r] (A2) 

where [cQg is the mxm matrix in the upper left-hand corner of [pis* 

The last term on the right is the deflection of the supports calculated 
by applying the (equal and opposite) reaction forces on the fuselage. 

From equation (A2), 

-([j] s + DS s )w - -MW - + eW + e[rl (as) 

where now jV] = [j]g + [g] s is a symmetric mxm matrix. The inverse 

of [t] must be confuted, but it is ordinarily a low-order matrix. From 
equation (A5), the following equation is obtained: 



NACA TN 36kO k-9 

-W - W^MgH + + e [ t] - 1 [r] (a4) 

Inasmuch as the wing is in equilibrium under the loads [p] and [jr], 

the resultant vertical force and pitching moment due to these loads must 
vanish; this condition yields the following two equations : 

[I]'[e] + [i]'[p] - 0 («) 

M'M + H'H-o w 


where primes denote the transpose of a matrix. Combining equations (a 4) 
to (a6) gives the following two equations for the determination of £ 
and 0 : 

H 'H - [r] , [i]- 1 [5] s [p] + cM'fr]- 1 ^] + e[r]'[i]-i[r] (A7) 

W'H ■ DG'tfflsH ♦ ♦ •M’&Mfl w» 

or 

[t] ' f t] - 1 jr] 5 + [x] ’ [t] ■ - 1 jj] 9 - ( [l] ' - [r] ' [t] - 1 [a] g ) [p] ( A9 ) 

[r]'[T]^[r]e - ([!]' - [r]'W-[oj s )[p] ( «o) 

These equations must now be solved for £ and 0. The following (scalar) 
numcers have to be computed: 
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(W'Hl'Mfl'HIrl - ([r]'[.ra 

• - iwm* 

• (WHIr])/. 

■ (WH-'wy* 


(All) 


Solving for £ and 0 gives 




• ■ 4]' - [r] 1 H1i)H - 

lj - [t] [tJ [g s j P_ 


(A12) 

(A13) 


With the help of the matrix |t] b _1 ^ which is [t] - 1 bordered by zeros 
■to make up an MXM matrix') , these equations may be written : 


S - (-[!]■ - >[«]’) (W - Hb’U)H 

e - (r[ { ]' - nH*) (H - 


(A14) 

(A15) 


where is the unit matrix. Now that £ and 8 

may be found from equation (a4) and is as follows: 


are known, 



r 
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Introducing these results into equation (Al) gives 

M - (H - HsHb' 1 ) {[o]s + (“MH' - pHW - 
pHH' + ?&] [e]') (CO - Hb-HsJH ' a1 ?> 

or 

M - [ K J s [ p ] ^ 

The square symmetric matrix jj<TJg whose definition is obvious from equa- 
tions (A17) and (A18) is the required modified set of influence coeffi- 
cients for the wing. 


Results in the Antisymmetric Case 


The modified matrix of influence coefficients in the antisymmetric 
case is simpler because only one rigid-body motion is involved, namely 
a roll. In this case. 


*Ja - (CO - HaBb- 1 ) 



+ 



(Ai9) 


where now [t[] is the column of tj's at the first m reference points 
and f Til is the column of t] ' s at all M reference points . Also 


and 


M- [ J ]a + [o]a (A 20 ) 

- WBBO 


1 

a 


(A21) 
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Exceptional Support Conditions 

In general the jjrj matrix in equation (A3) will he nonsingular, 

hut in certain special cases it will he singular and some modification 
of the method will he necessary. These cases are listed as follows: 

(1) The [g] matrix will he singular if there are more than three 

supports along any one chord in the parabolic theory (the number is two 
in the linear theory) because the influence coefficients for four or more 
points along a chord which deforms into a parabola are linearly dependent. 

If the supports are fixed ^so [ jJ = 0^, then £tJ will he singular. 

There is no remedy for this situation other than to reduce the number of 
supports - and this reduction does not result in any loss of generality 
within the framework of parabolic theory. 

(2) In the symmetric case, if the origin is to he a fixed support, 
then exclude the origin from the list of reference points. Since the 
origin is fixed, £ will he zero and one equation must he deleted from 
the set of equations used to determine all the unknowns; this will he 

equation (A5). The effect on the result, namely [k]c, in equation (Al8), 

is to make 



(3) In the symmetric case of the linear theory, if any point along 
the center line is to he fixed, then exclude this point from the list of 

reference points. For purposes of forming the and matrices, 

choose a new y-axis passing through this point and set £ = 0. Again the 
effect on the result is to make 


a =3 = 0 7= 

(4) In the antisymmetric case, if the points (x=0,y=±e) are to be 
fixed, then exclude these points from the list of reference points and 
set a = 0 in equation (A19)- 


Miscellaneous Remarks 

In cases where some of the points of attachment cannot conveniently 
be made to fall -on a station, influence coefficients for these points can 
he obtained by interpolation from surrounding reference points. 
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The effect of a couple applied at a point of attachment can be 
roughly approximated by replacing the couple by equal and opposite forces 
at nearby reference points . 

Clearly, a set of influence coefficients for the whole fuselage 
could be modified in exactly the same way as in the previous section in 
order to take account of the presence of the wing. From there it would 
be only a step to obtain an overall set of influence coefficients for 
the wing-fuselage combination, the connecting link being furnished by the 
expressions for the reactions at the points of attachment. Since the 
required matrix manipulations are quite elementary, and the present paper 
is concerned primarily with the wing, the details will not be given herein. 

A set of influence coefficients obtained experimentally could be 
modified by the method given in this appendix to obtain a modified set 
appropriate to same other system of supports than that used during the 
test. The only provision is that the support system used in the test be 
such that the reactions at the supports would vanish whenever the struc- 
ture is subjected to a self -equilibrated set of loads. Fixing the struc- 
ture at three points is one possibility. 
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APPENDIX B 


EXAMPLE DERIVATION OF A MATRIX EQUATION 

Sn g 

In this appendix the matrix form of the quantity is derived. 

° c Pk:,n 

This quantity is the term contributed by a spar to equations (28). The 
derivation will be made for the symmetric case . 

From equation (22) 

= ^(h,o*a,o S + fs.iVi 2 + • • • + (B1) 


where 


^s,n " x s,n-l " ^ x s,n + x s,n+l 


(B2) 


x s,n = *fio,n + x s,nfPl,n + 


(B3) 


The boundary conditions are (from eq. (25)) 


* 0,0 = *** 1,0 “ ° * 0,-1 “ * 0,1 * 1,-1 " * 1,1 


The first few of equations (B3) read 


* 2,-1 = * 2,1 


(b4) 


*S,-1 “ <1*0,1 + x s, -1*1,1 + x s,-l 2 <P2,]J 


X s,0 = x s,0 *2,0 


X s,l = *0,1 + X s,l*l,l + X s,l *2,1 
X s,2 = ^0,2 + *6, 2*1,2 + x b,2 2 $2,2 


> 


(B5) 
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The first few of equations (B2) read 

^s,0 = *s,-l “ ^SjO + *8,1 



= X - 

- 2X n 

+ 

X ^ 

S,1 

s,0 

S,1 


s,2 


= X , 

- 2X ^ 

+ 

X 

s,2 

8,1 

s ,2 


S,5 


(b6) 


The first two equations of ~^ s read (with the use of eqs. (B5) aad 

^Po.n 

(B6) ) 


^7 = ^3 ( 2p 6,0+s,0 “ 2Ps,l+s,l + Ps,2V2) 


Scp 0 
Sn 


Sep, 


' 0,2 


= ^(Ps, 1+8,1 - 2 Ps,2+S,2 + Ps, 3+6,3) 


(BT) 


gj-j- 

The first two equations of - — — read. 


^l,n 


- -^[( X 6,-l + ^l) ^S, 0+8,0 “ 2X S,lPs, 1+8,1 + X B,lPs,2+B,2j 


^XX / 

dep^g = ^( x s,2Ps,l+s,l - , 2^*8 , 2+s , 2 + x s, 2^s, 3+8,3) 

an. 


The first two equations of 


Sep 


read 


2,n 


- * s , 0 %, 1 *° a ) 

an, 

&P; 


*7 = ^3 [( x s,-l 2 + X s,1 2 )3s,0+b,0 - 2x s,1 2 ?b,*,1 + x s,1 2 Ps,2+s,2 ] 


(B8) 


(B9) 
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Freon equations (B7) , the first two rows of the matrix 
written as follows: 


an s 

*?0,n 


may he 


fans 


rH 

e? 

so 


3n s 

_ 1 

* 

“ e 3 

• 

« 



-2 1 
1-2 1 


Ps,o 


s,l 


*s,0 

*B,1 


(BIO) 


Similarly, from equations (B8), 


dn s 

^,1 


* 8,-1 + * 8,1 

- 2 * 8,1 * 6,1 

Ps,0 



* 8,0 

Sn a 

** 1,2 

X 

' e 3 


* s,2 - 2 * 8,2 * 8,2 

* S,1 



* 8,1 

, 







• 


and, from equations (B9), 


SlTjj 


- 2 * 8 , 0 2 

* s , 0 2 



^ 8,0 



^ 2,0 





* 8,0 

Sn B 

1 

X ,2 . x 2 

-2*8, l 2 

X 2 

* 8,1 


^ 8,1 


* 8,1 

^> 2,1 

m 7 

* 8,-1 + * 8,1 



• 



* 8 , 2 2 

Py 2 v 2 

- 2 * 8,2 * 8,2 


. . . 


• 

• 




• • 


J 


• 


( Bll ) 


( B12 ) 


Obviously, if more rows were written the 1-2 1 pattern would he 

followed until the last few rows, where the vanishing of p for 

8;I1 

n > N changes the pattern. The first few rows of j~t s may he 
written 
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Vo 


2 


Vi 


x s,-l + x s,l 



\,1 

Vi 


-2 1 


90,2 

+ 

' 2x b,1 

*6,2 


'Pi, 2 

• 


1-2 1 




x s,l 

-2*8, 2' ^,3 


* 

• 


■ • • 


* 



. . . 


• 


-^s.O 2 x s,-l 2 + x s,l 2 
x s,0 2 -2x s,l 2 


x s,2 2 

_2x Sj2 2 x s,3 2 


«P. 


2,0 



(B13) 


The three rectangular matrices appearing in equation (B13) are the 
transposes of the three rectangular matrices appearing in equations (BIO), 
(Bll) , and (B12) . The form of these matrices in their lower right-hand 
corners can he obtained in a way similar to that which was just us^d to 
obtain the form in the upper left-hand comer. Equations (BIO), (Bll), 
and (B12) can be written compactly as follows: 



These equations can be further combined to give 


Sq> 0 


WPJW MMH' 

WWW 

'Po 

to. 

_ 1 
e 3 

mm mh&t 

WWW 

*1 

aiie 

Sqig 


MP-1W MHW 

MMM 



A comparison of equation (B17) with equation (30) will show how the sth 
spar contributes to the matrices defined in equations (31). Simi- 

lar derivations were made for the covers and ribs but because the algebra 
is long and tedious only the results are given in the text. 
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(b) Antisymmetric. 

Figure 4.- Symmetric and antisymmetric support conditions considered. 
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